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ABSTRACT 



• We present the first high-precision photometry of the transiting extrasolar planetary system WASP-7, obtained using telescope defo- 

• ' cussing techniques and reaching a scatter of 0.68 mmag per point. We find that the transit depth is greater and that the host star is more 

(N ■ evolved than previously thought. The planet has a significantly larger radius (1.330 ± 0.093 Rjup versus 0.915^q 55^q Rjup) and much 

' — ^ . lower density (0.41 ± O.lOpjup versus 1.26^55 2[ pjup) and surface gravity (13.4 + 2.6 m s"^ versus 26.4;^^ jj m s"^) than previous mea- 

' surements showed. Based on the revised properties it is no longer an outlier in planetary mass-radius and period-gravity diagrams. 

We also obtain a more precise transit ephemeris for the WASP-7 system. 

Key words, stars: planetary systems — stars: individual: WASP-7 

■><■ 

. ! 1 . Introduction terisation of WASP-7 therefore relied upon the photometry ob- 
tained by the WASP-South telescope tPollacco et al. 2006) . 

The transiting extrasolar planet (TEP) WASP-7 b was discovered 
by the WASP consortium (Hellier et al. 2009, hereafter H09) 

through the detection of transits in front of its F5 V pai-ent star. It The relatively large scatter in the discovery data meant that 

is a challenging target for acquiring high-precision transit pho- the transit shape was poorly delineated. Because of this, the 

tometry, due to the brightness of the parent stai- (V = 9.5), the analysis by H09 included an additional constraint in the form 

paucity of good nearby comparison stars, the transit duration of a main seq uence mass-radius relation for the host star (e.g. 

(3.8 hr), and the orbital period (4.95 d) which is both compar- lAnderson et al. 2010) . The radius, surface gravity and density of 

atively long and close to an integer number of days. The charac- the planet resulting from their analysis are R^, = 0.91!° °^[; Rjup, 

gb = 26.4;^^yms"^ and pb = l-26;^|j 2j pjup, respectively. These 

* Based on data collected by MiNDSTEp with the Danish 1 .54m tele- values placed WASP-7 b in an outlier position in the mass-radius 

scope at the ESO La Silla Observatory. diagram of TEPs, having one of the largest densities within the 

** Royal Society University Research Fellow main planet population (masses $2 Mjup). This was interpreted 
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by H09 as evidence that WASP-7 b has a massive heavy-element 
core. 

In this work we present the first follow-up photometric ob- 
servations obtained for WASP-7. The high precision of our ob- 
servations (0.68 mmag scatter) allows us to obtain the physical 
properties of the transiting system without needing to impose 
any constraints on the parameters of the parent star We find a 
substantially larger radius, and therefore a lower density and sur- 
face gravity. We also greatly improve the orbital ephemeris for 
the system, so transit midpoints in the 2011 observing season 
can be predicted to within 45 s instead of 27 min. 

2. Observations and data reduction 

WASP-7 was observed from 00:14 to 07:18 UT on the date 
2010/09/07. We used the 1.54m Danish telescope at ESO La 
Silla, equipped with the DFOSC focal-reducing CCD imager 
A total of 217 integrations were obtained through a Gunn / fil- 
ter (ESO filter #425) and with an exposure time of 60s. The 
CCD was windowed down in order to reduce readout time, re- 
sulting in a dead time of 50 s between individual observations. 
The telescope was defocussed to a point spread function diam- 
eter of 38 pixels (15") in order to average out flat-fielding noise 
and avoid saturating the CCD. 

The data were reduced using a pipeli ne which i mplements 
the DAOPHOT aperture photometry routine (,Stetsonll987 ). The im- 
ages were debiassed and flat-fielded, and slight pointing vari- 
ations were detected and accounted for by cross-correlating 
against a reference image. A differential-magnitude light curve 
was constructed by simultaneously fitting a straight line to the 
out-of-transit observations plus weights to the reference stars 
used to create the ensemble comparison star. We find that the 
sizes of the software apertures and the choice of comparison 
stars has very little effect on the shape of the observed transit. For 
our final photometry we adopt the light curve which displays the 
lowest scatter with respect to a fitted transit model (see below). 
The root-mean-square (rms) of the scatter in this light curve is 
0.68 mmag per point. For furt her details on the approach used to 
obtain and reduce the data see lSouthworth et a"hC2009., .201ft) . 

3. Light curve analysis 

The light curve of WASP-7 was modelled using the jk- 
TEBOP codfl which in turn is based on the ebop program 
(Popper & Etzel 1981; Nelson & Davis 1972). The main param- 
eters of the fit were the sum of the fractional radii of the star 
and planet, ta -H rt = ^hlEh. where a is the semimajor axis, their 
ratio, k - J2 ™d the orbital inclination, After a preliminary 
fit was obtained, the measurement errors were rescaled to give 
a reduced x' of 1 .0. In order to provide an accurate measure 
of the orbital period we included the time of mid-trans i t quote d 
by H09, using the approach given bv lSouthworth et alJ (l2007ah . 
This results in an orbital ephemeris of: 

To = BJD(TDB) 2455 446.63493(30) + 4.9546416(35) x £ 

where E is the number of orbital cycles after the reference epoch 
and the bracketed quantities denote the uncertainty in the fi- 
nal digit of the preceding number This ephemeris is on the 
BJD(TDB) timescale (Eastman et al. 2010). The timestamps in 
our data were manually verified to be correct to within +2 s. 

' JKTEBOP is written in fortran77 and the source code is available at 
|http : //www ■ astro . keele . ac . uk:/~ jkt/codes/ j ktebop . html , 



Table 2. Final parameters of the fits to the light curve of WASP- 
7, compared to the results from H09. 



This work H09 



'■a + 'b 


0.1207 


± 0.0068 


0.1001 
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0.0956 


±0.0016 


0.0761 ± 0.0008 
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±0.93 


89 6+° '* 




0.1102 


±0.0061 


0.0930 


rh 
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± 0.00070 


0.00708 
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Fig. 1. Phased light curve of WASP-7 compared to the best fit 
found using jktebop. The residuals of the fits are plotted at the 
base of the figure, offset from zero. 



Limb darkening (LD) was included and its uncertainty ac- 
counted for using five different functional laws, and solutions 
were obtained for the cases when both coefficients were fixed to 
theoretical values, the linear coefficient was fitted for, and both 
coefficients were fitted for Theoretical values were calculated 
for the measured effective temperature and surface gravity of the 
star by linearly interpolating within the tables of ' Van Hammd 
(Il993l) andlClaret (2000). The results for each of these modelling 
runs are given in Table[T] We find that the results for different 
LD laws are generally in good agreement, although the amount 
of LD is lower than theoretically predicted. 

Uncertainties in the fitted parameters were obtained using 
Mont e Carlo and residual- permutation simulations ( SouthworthI 
200i Llenkins etalJl2002h . We are able to reject the solutions 
with all LD coefficients fixed to theoretical values, as they pro- 
vide a clearly worse fit to the data. The solutions where two LD 
coefficients are fitted are marginally the best ones, so for the fi- 
nal parameters we adopt their mean values from these four so- 
lutions. We adopt the parameter uncertainties from the residual 
permutation algorithm, which are roughly 1.5 times larger than 
the Monte Carlo alternatives. This implies that correlated noise 
is significant in our data, as expected from visual inspection of 
the plot of the best fit (Fig.[T]). The final parameter values and 
uncertainties are given in Table|2l Compared to those from H09, 
we find a solution with a lower orbital inclination and larger frac- 
tional radius for both the star and the planet. Our results were 
obtained using the methods of ISouthworthI (l2008h and are ho- 
mogeneous with the results in that work. 
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Table 1. Parameters of the jktebop best fits of the light curve of WASP-7, using different approaches to Umb darkening (LD). 



Linear LD law Quadratic LD law Square-root LD law Logarithmic LD law Cubic LD law 
All LD coefficients fixed: 



/-A + Tb 0.1255+0.0052 0.1180 + 0.0049 0.1194 + 0.0048 0.1184 + 0.0051 0.1223+0.0046 

k 0.09598 + 0.00089 0.09480 + 0.00075 0.09524 + 0.00075 0.09475+0.00082 0.09633+0.00060 

iO 86.61+0.52 87.43 + 0.59 87.23 + 0.56 87.39 + 0.63 86.80 + 0.47 

UA 0.43 fixed 0.20 fixed 0.00 fixed 0.55 fixed 0.20 fixed 

Va 0.30 fixed 0.60 fixed 0.25 fixed 0.15 fixed 

rA 0.1145 + 0.0047 0.1078 + 0.0044 0.1090 + 0.0044 0.1081 +0.0046 0.1116 + 0.0041 

T-b 0.01099 + 0.00053 0.01022 + 0.00048 0.01038 + 0.00048 0.01025 + 0.00051 0.01075 + 0.00045 

tr (mmag) 0.7136 0.6977 0.6946 0.7015 0.6828 

Fitting for the linear LD coefficient and perturbing the nonlinear LD coefficient: 

/-A + fb 0.1224 + 0.0043 0.1192 + 0.0047 0.1227 + 0.0046 0.1221 +0.0047 0.1200 + 0.0044 

k 0.09675+0.00069 0.09587 + 0.00072 0.09646 + 0.00065 0.09628+0.00069 0.09623+0.00071 

86.77 + 0.44 87.15 + 0.55 86.77 + 0.48 86.83 + 0.49 87.05+0.49 

UA 0.245 + 0.047 0.059 + 0.064 -0.134 + 0.059 0.387 + 0.064 0.190 + 0.051 

Va 0.30 perturbed 0.60 perturbed 0.25 perturbed 0.15 perturbed 

rA 0.1116 + 0.0039 0.1088 + 0.0043 0.1119 + 0.0042 0.1113+0.0043 0.1095+0.0039 

rfc 0.01080 + 0.00043 0.01043 + 0.00047 0.01079 + 0.00045 0.01072 + 0.00046 0.01054 + 0.00042 

cr (mmag) 0.6841 0.6829 0.6826 0.6825 0.6827 

Fitting for both LD coefficients: 

TA + rb 0.1216 + 0.0050 0.1199 + 0.0049 0.1221 +0.0054 0.1196 + 0.0050 

k 0.09583 + 0.00115 0.09537 + 0.00145 0.09575 + 0.00148 0.09527+0.00151 

iC) 86.91 +0.56 87.15 + 0.60 86.89 + 0.63 87.18 + 0.63 

Ma -0.021 +0.256 -0.866 + 0.979 0.578 + 0.391 0.066 + 0.131 

Va 0.42 + 0.41 1.83 + 1.60 0.56 + 0.60 0.54 + 0.45 

Va 0.1110 + 0.0045 0.1095 + 0.0045 0.1114 + 0.0049 0.1092 + 0.0045 

rb 0.01064 + 0.00051 0.01044 + 0.00051 0.01067 + 0.00055 0.01040 + 0.00052 

cr (mmag) 0.6826 0.6818 0.6820 0.6817 



Table 3. Derived physical properties of WASP-7 using five different sets of theoretical stellar models. 





{Claret models) 


(Y^ models) 


(Teramo models) 


(ra55 models) 


{DSEP models) 


2fb (kms-') 


135.6 + 2.1 


135.6+1.9 


135.1+2.2 


135.4+1.9 


134.4 + 2.2 


Ma (Mo) 


1.285 + 0.059 


1.287 + 0.054 


1.273 + 0.061 


1.280 + 0.052 


1.254 + 0.061 


i?A(Ro) 


1.436 + 0.092 


1.436 + 0.092 


1.431 +0.092 


1.434 + 0.090 


1.424 + 0.092 


loggA (Cgs) 


4.233 + 0.047 


4.233 + 0.047 


4.231+0.047 


4.232 + 0.047 


4.229 + 0.047 


Mb (Mj„p) 


0.96 + 0.13 


0.96 + 0.13 


0.96 + 0.13 


0.96 + 0.13 


0.95+0.13 


(Rjup) 


1.333 + 0.093 


1.334 + 0.093 


1.329 + 0.093 


1.331 +0.092 


1.322 + 0.092 


Pb (Pjup) 


0.407 + 0.101 


0.406 + 0.101 


0.408 + 0.101 


0.407 + 0.101 


0.410 + 0.102 





0.069 + 0.011 


0.069 + 0.010 


0.070 + 0.011 


0.070 + 0.011 


0.070 + 0.011 


a(AU) 


0.06185 + 0.00095 


0.06188 + 0.00086 


0.06166 + 0.00098 


0.06177 + 0.00084 


0.06134 + 0.00100 


Age (Gyr) 


2 5+"-^ 


2 5+"- 


9 -1+0.6 




9 7+0.8 



4. The physical properties of WASP-7 

The light curve parameters alone do not allow the physical prop- 
erties of the WASP-7 system to be calculated. Additional ob- 
served properties are available from H09: the star's orbital ve- 
locity amplitude (Ka = 97 + 13 ms"'), effective temperature 
(Feff = 6400 + 100 K), surface gravity (logg = 4.3 + 0.2) and 
metallicity ([§] - 0.0 + 0.1). Of these, the surface gravity is 
quite imprecise so is only useful as a consistency check. 

In order to calculate th e physical properti es of WASP-7 we 
use the same approach as in ' S out h worthi (l2009l) . Firstly, we adopt 
tabulated predictions from theoretical models of low-mass stars. 
Secondly, we guess an initial value of the velocity amplitude of 
the planet, K^,. We then use the known ta, r\,, i, Ka and orbital 
period to calculate the physical properties of the system. The ob- 
served Tgff and calculated Ra are then compared to the model- 
predicted values for a star of the calculated mass (Ma), and K\, 
is adjusted until the best match is found. This is performed for 
a grid of ages covering 0. 1 to 20 Gyr in 0. 1 Gyr steps, by which 



the overall best physical properties and age are found. The uncer- 
tainty in each of the input parameters is propagated by repeating 
this procedure with the parameter perturbed, resulting in a de- 
tailed eiTor budget for every output parameter. 

The above procedure is performed five times, using in- 
depende nt sets of th eoretic al model predict i ons. These are: 



Claret dClaretl l2004 dDemarque et al.' '2004), Teramo 



jPietriiiferni et al.l |2004. VRSS tVandenBerg et al.. .200 6) and 
DSEP (iDotter et al.ll2008h . The final value for each of the phys- 
ical properties is taken to be the unweighted mean of the values 
found using the five different model sets. The random error is 
taken to be the largest of the individal errors found by the pertur- 
bation analysis. The systematic error is the standard deviation of 
the values from the five different model sets. The results for each 
model set are given in Table[3] The final physical properties of 
WASP-7 are collected in Table|4] and compared to those found 
by H09. In this table the quantity T^'^ is the equilibrium tem- 
perature of the planet excluding the energy redistribution factor 
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Table 4. Final physical properties of WASP-7 compared to those 
found by H09. The second eiTorbars represent systematic errors. 





Thi*; work" Tfirial^ 


H09 


Ma (Mr,) 


1.276 ± 0.061 +0.022 


1.28+""^ 




1 + 092 + 008 


1 236+°'''^''' 


6 A V,^&^^ 


4 2M + 047 + 003 


4 363+'' '''° 


Pa (Po) 


0.434 ± 0.074 


Mb (Mj„p) 


0.96 ± 0.13 ±0.01 


96+" '^ 


(Rjup) 


1.330 ± 0.093 ±0.008 


n qi ^+" 046 


gb (ms"') 


13.4 ± 2.6 


26 4+" 


Pb (Pjup) 


0.41 ±0.10 




T^, (K) 


1487 ± 48 


1379!^^ 





0.070 ± 0.011 ±0.000 




a (AU) 


0.0617 ±0.0010 ±0.0004 


0.0618!!J;™^ 


Age (Gyr) 


7 4 +0.8 +0.3 
^■^-10-04 




1 .0 

Planet mass 



Fig. 2. Plot of the masses and radii of the known TEPs. The 
blue symbo l s deno te values from the homogeneous analysis of 
ISouthworflJ ( 12010 ) and the red symbols results for the other 
known TEPs. WASP-7 is shown in black with an open circle 
(H09) and a filled circle (this work). Grey dotted lines show 
where density is 1.0, 0.5 and 0.25 pjup. 



(ISouthworthll2009h and is the ISafronovl(ll972h number. We 
find that the star is slightly evolved, allowing a reasonable esti- 
mate of its age (2.4 + 1.0 Gyr). 

5. Summary and conclusions 

We have observed a transit of the WASP-7 system using defo- 
cussed photometry techniques, resulting in a light curve with an 
rms scatter of only 0.68 mmag. These data were modelled with 
the JKTEBOP code, and the results used to determine the physical 
properties of the planet and its host star. The higher quality of 
our photometry has allowed us to measure the system proper- 
ties without constraining the host star to follow a main-sequence 
mass-radius relation. Compared to H09, we find both a larger ra- 
dius for the host star and a larger ratio of the radii of the planet 
and star. This results in a sizeable increase in the planetary ra- 
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Fig. 3. Plot of the orbital periods and surface gravities of the 
known TEPs. The symbols are the same as for Fig.|2] 



dius, from 0.915!j}°46Rj^p i 339 + 0.093 Rjup, which in turn 
means a lower surface gravity and mean density. 

We have collected literature measurements of the physical 
properties of the 96 known TEPs as of 2010/11/18. For 30 of 
these objects we used the results of the homogeneous analysis 
performed by Southworth (2010). In Fig.|2]we plot their masses 
and radii, plus the values from WASP-7 found in this work (filled 
circle) and from H09 (open circle). The H09 results for WASP- 
7 b indicated that it was one of the densest known TEPs be- 
low ~2 Mjup, along with the more recently-discovered planets 
CoRoT-13 b (Cabrera et al. 2010) and HAT-P-15 b (Kova cs et al.1 
201(J). The outlier status has now been lifted: our results place 
WASP-7 b in a well-populated part of the mass-radius diagram 
and demonstrate that high-quality data is required to obtain reli- 
able measurements of the properties of TEPs. 

Fig.|3] shows the surface gravities of the known TEPs as a 
function of their orbital periods. The existing planet population 
shows an inverse correlation between period and surface gravity 
(Southworth et al. 2007b), at least for the domina nt population 
with periods $10 d and masses $3Mjup (see also iFressin et al.l 
I2009D . The revised properties of WASP-7 b move it from outlier 
status to within the sprawl of parameter space occupied by the 
general planet population. 

The theoretica l mo d els o f irradiated gas giant plan- 
ets Bodenheimer et al.1 (l2003h . iFortnev et 5] (§007) and 
iBaraffe et al., (i2008]) predict radii of no more than 1.16 Rjup for 
a l.OMjup planet with a range of chemical compositions and 
core masses, and without an arbitrary additional heating source. 
Our upward revision of the radius of WASP-7 b means that it no 
longer matches the predictions of these models at the 2cr level. 
These conclusions could be strengthened by the provision of 
more radial velocity measurements as well as the acquisition of 
further photometric observations. 
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